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We present explicit expressions for the optical scalars and the deflection angle in terms of the 
energy-momentum tensor components of matter distributions. Our work generalizes standard ref- 
erences in the literature where normally stringent assumptions are made on the sources. 



I. INTRODUCTION 

Gravitational lensing has became a significant tool 
to make progress in our knowledge on the matter con- 
tent of our Universe. In particular, there is a large 
number of works that use gravitational lensing tech- 
niques in order to know how much mass are in galaxies 
or clusters of galaxies. One of the most exiting results 
was to reaffirm the need for some kind of dark matter, 
that appears to interact with the barionic matter only 
through gravitation. 

The question in which there is yet not general agree- 
ment is on the very nature of this dark matter. The 
most common conception is that it is based on colli- 
sionless particles WeiOjB] , and where the pressures are 
negligible. However in the context of cosmological 
studies, one often recurs to models of dark matter in 
terms of scalar fields. There is also the possibility that 
dark matter were described in terms of spinor fields. 

One method to study the nature of dark matter con- 
sists in observing the deformation of images of galax- 
ies behind a matter distribution that is the source of 
a gravitational lens. 

In many astrophysical situations, the gravitational 
effects on light rays is weak, and the source and ob- 
server are far away from the lens, therefore they are 
studied under the formalism of weak and thin gravita- 
tional lenses. The basic and familiar variables in this 
discussion are shown in figure [I] 

In this framework the lens equation reads 

I3 a ^6 a -^a a . (1) 

The differential of this equation can be written as 
5(3 a = A\ S9 b , (2) 
where the matric A a b is in turn expressed by 

A a _ ( 1 - « - 71 "72 V ,„v 

Ab ~{ -72 1-K + 7 J' [) 

where the optical scalars k, 71 and 72, are known 
as convergence k and shear components {71,72}, and 
have the information of distortion of the image of the 
source due to the lens effects. 




FIG. 1. This graph shows the basic and familiar angular 
variables in terms of a simple flat background geometry. 
The letter s denote sources, the letter I denotes lens and 
the observer is assumed to be situated at the apex of the 
rays. 

It is somehow striking that in most astronomical 
works on weak gravitational lensing, it is assumed that 
the lens scalars and deflection angle, can be obtained 
from a Newtonian-like potential function. These ex- 
pressions although are easy to use, have some limita- 
tions: 

• They neglect more general distribution of 
energy- momentum tensor T a b, in particular they 
only take into account the timelike component 
of this tensor. In this way they severely restrict 
the possible candidates to dark matter that can 
be studied with these expressions. 

• They are not expressed in terms of gauge invari- 
ant quantities. 

• Since these expressions are written in terms of a 
potential function, it is not easily seen how dif- 
ferent components of T a b contribute in the gen- 
eration of these images. 

Moreover, most of them assume from the beginning 
that thin lens is a good approximation. 

In other cases in which the thin lens approximation 
is not used[BBV10 , the results are presented in a way 
in which gauge invariance is not obvious. 



2 



Here we extend the work appearing in standard 
references on gravitational lensing [SEF921 ISSE941 
IWammbj IBarlOj and present new expressions that do 
not suffer from the limitations mentioned above. In 
particular we present gauge invariant expressions for 
the optical scalars and deflection angle for some gen- 
eral class of matter distributions. 



II. INTEGRATED EXPANSION AND SHEAR 

A. General equations: The geodesic deviation 
equation 

Let us consider the general case of a null geodesic 
starting from the position p s (source) and ending at 
p (observer). Let us characterize the tangent vector 
as I — t^; so that 



0: 



(4) 



that is, A is an afhne parameter. 

We can now consider also a continuous set of nearby 
null geodesies. This congruence of null geodesies can 
be constructed in the following way. Let S be a two 
dimensional spacclike surface (the source image) such 
that the null vector £ is orthogonal to S. Next we can 
generalize £ to be a vector field in the vicinity of the 
initial geodesic in the following way: let the function 
u be defined so that it is constant along the congru- 
ence of null geodesies emanating orthogonally to S 
and reaching the observing point p . Then, without 
loss of generality we can assume that 

4 = V a u; (5) 

which implies that the congruence has zero twist. 

We can complete to a set of null tetrad, so that m a 
and m a are tangent to S. Then a deviation vector at 
the source image can be expressed by 

(6) 



In order to propagate this deviation vector along the 
null congruence one requires, that its Lie derivate van- 
ishes along the congruence; that is 

- 0. (7) 

Using the GHP notation[GHP73] it can be shown 
that the previous equation can be written as 

D(c)+Cp + C* = 0; (8) 

where P is the well behaved derivation of type {1, 1} 
in the direction of £ (the null geodesic vector of the 
congruence). 
Defining X by 



X 



one can prove that ^ satisfy 

£(£(X)) = -QX 
where Q is given by 

' $_oo *o 
*o $00 



Q = 



(9) 
(10) 

(11) 



with 



and 



*o = C abcd £ a m b £ c m d . 



(12) 



(13) 



Therefore, this form of the equation only involves cur- 
vature quantities. 

Although the last equation can be integrated nu- 
merically without problems; it is sometimes conve- 
nient to have at hand some method for approximated 
solutions. So, next we present an approximation 
scheme that it can be applied to any order one wishes 
to obtain; although we will concentrate on the linear 
approximation since in weak lens studies it is con- 
sistent to consider linear effects of the curvature on 
geodesic deviations. 

Let us first transform to a first order differential 
equation. Defining V to be 



and 

one obtains 

*(X) 

with 



V = 



X = 



dX 
~dX~'' 



dX 



A = 



V 

-QX 



I 
-Q 



AX; 



(14) 



(15) 



(16) 



(17) 



Equation ( 16 1 can be reexpressed in integral form, 
which gives 



X(A) = X +/ j4(A') X(A') dX'. 

An 



One can define the sequence 



Xi(A) =X + f A(X')X dX', 



(18) 



(19) 



X 2 (A) = X +/ AWXxWdX'; (20) 
and so on. 

Assuming that Q is in some sense small, one expects 
that this sequence will converge and therefore provide 
for the solution. 

The complete linear iteration is given by 
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X 3 (A) = 



1 - fL fx, \ Q" dX " dX ' ( A - A o)! - ft Ix' ( A " - W d\" dX' 



' Ao ^Ao 



I-/£(A'-Ao)g'dA' 



(21) 



Now in order to integrate the geodesic deviation equa- 
tion, we must choose the correct initial conditions. In 
the case of light rays belonging to the past null cone 
of the observer and intersecting S at the source, this 
initial conditions are X = and V ^ 0; since one 
can think the beam, starts backwards in time from 
the observer position, and so initially has vanishing 
departure, but with nonzero expansion and shear. 
Therefore in the linear approximation one has 

*(A)=((A-A )I 

" A \ (22) 

(A-A')(A'-Ao)Q'dA' V(A ); 1 ' 



But in the case of the presence of a gravitational 
lens, if an observer sees an image of "size" X Q , which 
means X = d s V Q (since actually what is observed is 
V = V(A ) ), then it should be produced by a source 
of size X s = X(X S ), as described by equation (22 1. 

We set from now on Ao = and A s = d s , then 



eq.(22) reduces to 



X' 



l s JQ 



X'{d s - X')Q' dX' )x o . (25) 



B. Optical scalars in terms of curvature 



and 



V(A) =y- j x ( A ' - W dX'j V(A ). (23) 

In these integrations Ao indicates the position at the 
observer and from now on, A s will indicate the position 
at the source. 

We observe from the first expression, that if the 
metric were flat (Q = 0), in order to get a deviation 
vector constructed from X\ , defined as X evaluated at 
A s = A + d s , one must choose as initial condition 



V(A ) = 1 , , X(X S = X +d s ) = ^-Xl 
[A s — ao) a s 



In order to compare with the standard representa- 
tion of the lens scalar we note that the original de- 
viation vector in the source will be given by eq.(25) 
i.e. 




(26) 



if we make the following decomposition into real and 
imaginary part, 



Co — SoR 4 
*0 = *0fl 



ksi, 



(24) we obtain from eq.(26) that 



(27) 
(28) 
(29) 



W = I 1 - ^ Ko - <n) dx\ w - ^-^^dX'^, 



(30) 



Note also that in principle the integration must be 
made through the actual geodesic followed by a pho- 
ton in its path from the source to observer. How- 
ever the last expressions are valid only in the limit 
where the linear approximation is valid. If one con- 
siders a linear perturbation from flat spacetime, then 
the curvature components <I>oo and would be al- 
ready of linear order. Then, in the context of weak 
gravitational lensing, it is consistent to consider a 
null geodesic in flat spacetime; since the actual null 
geodesic can be thought as a null geodesic in flat 
spacetime plus some corrections of higher orders. 



Now, in order to compare with the usual expressions 
for the lens scalars k, 71 and 72, let us recall that they 
are defined via the relation eq.Q; but since it is a 
linear relation, one can relate the deviation vectors by 
the same matrix, namely 



(31) 



where {c*, are the spatial vector associated with 
{q s , <; } respectively. 



Therefore, by replacing into eq.(31), we obtain 

SsR = (1 - K - 7i)^ ofl - 72 W, (32) 

w = + (1 - k + 71 Ko/; (33) 
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which by comparing with eq. ( 30 1 , implies that 



1 



A'Of.-A'KodV, 



72 



<-s Jo 



X'ids-X'MjdA'. 



(34) 
(35) 
(36) 



Let us emphasize that these expressions for the weak 
lens quantities are explicitly gauge invariant, since 
they are given in terms of the curvature components. 
This is in contrast to the usual treatment of weak 
lensing found in the literature, which use for example 
equation (2.17) of reference |SEF92j as the source for 
the calculation of the lens scalars. 

Note that these expressions are valid for any weak 
field gravitational lens on a perturbed at spacetime, 
without restriction on the size of the lens compared 
with the other distances. As a final comment to this 
section it is important to remark that these expres- 
sions are valid for any weak gravitational lensing on 
a perturbed flat spacetime, without restriction on the 
size of the lens compared with the other distances. 
However, if we make use of the hypothesis of thin lens, 
these equations can be further simplified, as we will 
show below. 



III. THE THIN LENS APPROXIMATION 
A. The general case 

Now, we will consider the case of a lens whose size is 
small compared with the distances to the source and 
the observer. Let us consider a null geodesic coming 
from a source located at a distance d s from the ob- 
server, and at a distance di s from the lens, coming 
parallel to the y axis, but in the negative direction; 
we will use J to represent the impact parameter and 
$ to denote the angle of the trajectory as measured 
from the z axis, in the (z, x) plane. 

Then if we assume a thin lens, $oo and will be 
sharply peaked around A — di, where it is located and 
the expressions for the lens scalars are reduced to 



did 



on. 



where 



didi s ~ 
7i + *72 = — i— W , 



$00 = / $oodA, 
Jo 

*o = / *o^A, 
Jo 



(37) 



(38) 



(39) 



are the projected curvature scalars along the line of 
sight. 

We again emphasize that these expressions for the 
weak lens scalars are explicitly gauge invariants. 



B. The axially symmetric case 

1. The lens scalars in terms of projected Ricci and Weyl 
Scalars 

For axially symmetric lens , the projected curvature 
scalars are given by 

*oo(J) = / ' $oo(A')dA', (40) 



*o(J) = -e 2l UMJ)- 
where one can see that 



-2ii) 



* (A')dA'. 



(41) 



(42) 



The reason for the minus sign choice is that in 
many common astrophysical situations one would find 

MJ) > o. 



By replacing in eqs. (37) and (38), we obtain for 
the lens scalars 

disdi 



k = 



7i 



72 = 



d s 
dudi 



$oo(J), 
MJ) cos(2tf), 
MJ) sin(2i?), 



which implies that 



dudi j . 

7 = —i—MJ)- 



(43) 
(44) 
(45) 

(46) 



These equations can be compared to those of reference 
[FKN00 ; where they use different notation but similar 
content. 



2. Deflection angle in terras of projected Ricci and Weyl 
Scalars 



We wish now to express the deflection angle in terms 
of the curvature scalars. 

If we define the components of a 1 = (a 1 , a 2 ) as 



(a l ) = a(J)( Z j, X f); 



(47) 



then, it can be shown that the optical scalars can be 
written in terms of a( J) as 

1 di s di (da ajj y 
dJ ~~J 



7i 

72 



2 d s 
1 di s di 
2~dT 
1 di s di 
2~dT 



cos(2tf) 
sin(2tf) 



da a(J) 
dJ J~ 

da a{J) 
dJ J~ 



(48) 
(49) 
(50) 
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It is interesting to note that 

,„ q > • i r\ q\ didi s a{J) 
k - 71 cos(2i?) - 72 sin(2i?) = — — ; (51) 

from which, using eqs. ( 43||45 1, it is deduced that 

a(J) = J(*oo(J) + Mj)). (52) 

It is worthwhile to remark that this constitutes an 
equation for the bending angle expressed in terms of 
the gauge invariant curvature components in a very 
simple compact form. 

IV. DETAILED STUDY OF STATIONARY 
SPHERICALLY SYMMETRIC LENSES 

A stationary spherically symmetric spacetime can 
be expressed in terms of the standard line element 

ds 2 = a{r)dt 2 - b{r)dr 2 - r 2 (d9 2 + sin 2 6dip 2 ). (53) 

where it is convenient to define $(r) and m(r) from 

a(r) = e 2 ^ r \ (54) 

and 

1 



b(r) 



1 - 



2m(r) 



(55) 



For our purpose, it is more convenient to use a null 
coordinate system to describe the spherically symmet- 
ric geometry. Let us introduce then, a function 



u = t — r 



(56) 



where r* is chosen so that u is null. Then by inspec- 
tion of equation ( 53 1 one can see that 



dr* 

du = dt — dr = dt 

dr 



-dr: 



since then one has 
ds 2 = adu 2 + 2\fahdudr - r 2 (d0 2 + sin 2 8dip 2 ) 



(57) 



(58) 



It is natural to define the principal null direction lp 
from 



lp = du; 

which implies that the vector is 

We complete to a null tetrad with 

i) 



Up 



d_ 
du 



UK 



dr' 



(59) 



(60) 



(61) 



with the complex null vector 

V2P Q d 



mp 



r dC 



(62) 



in terms of the stereographic coordinate (. 
Therefore, one has 



ab 



and 



U = 



1 



26A 2 2 ' 



(63) 



(64) 



The more general distribution of energy- momentum 
compatible with spherical symmetry is given by 



T tt = ge 2 ^; 



T Pr 
rr ~ 2m ('') ^j 



Ptr z 



T w = P t r 2 sin(0) 2 ; 



(65) 



(66) 



(67) 



(68) 



where we have introduced the notion of radial compo- 
nent P r and tangential component Pt- 
The Einstein field equations 

G ab = -8nT ab , (69) 

in terms of the previous variables are 

dM 



dr 



dr 



= Aur g, 



47rr 3 R 



2M(r) 



(70) 



(71) 



(72) 



d^ d* \ 2M 

\dr 2 [ dr> J 1 r ' 

+ ^ 2 rf$ M dM 

dr r dr 

dM o 
-r—+M = 8nr 3 P t . 
dr 



One can show that the Ricci scalars defined with 
respect to the tetrad {£p,np,rhp,rhp} are given by 

47T 



nn = (g+Pr) 

a 



$n = n{g- P r + 2P t ), 
$ 22 = an(g + P r ) , 

A= ~(Q-P r -2P t ); 



(73) 

(74) 
(75) 
(76) 



G 



and the Weyl scalar, 



spherically symmetric physical fields, are given by 



47r , . m 

*2=f(Q-P r + Pt)--j- 
3 r J 



(77) 



Now, when one made a tetrad transformation be- 
tween the spherically symmetric tetrad and the tetrad 
adapted to the photon geodesic, one find that the 
function a{J) expressed in terms of the spherically 
symmetric null tetrad reads, 



At: did 



d s 
didi s 
d. 



d,. 



j2 

p+P r + —(P t -P r ) 



dy 



47r(p + P t - P r ) - 



3m 



(80) 



dy. 



These new expressions let us see explicitly the con- 
tributions of different components of the energy- 
momentum tensor on the optical scalars. Let us note 
that a couple of terms disappear in the isotropic case 
in which P r = P t . 



<t(J) = J 



3J 2 ~ 

^#2 



-2-(*n-7*oo) 



(78) 



in) 



dy. 



Note that in this case, the integration is on the coor- 
dinate y, instead of using arbitrary affine parameter. 
Also note that r 



VP- 



.'/- 



This constitutes an important explicit relation for 
the bending angle in terms of the curvature as seen in 
an spherically symmetric frame; which is the natural 
frame for the sources of the gravitational lens. 

In terms of the physical energy-momentum tensor 
it is obtained 



a{J) = J 



-d. 



3J2 



m(r) 4-7T 



g(r) 



(79) 



+4Tr{g{r) + P r (r))]dy. 



V. FINAL COMMENTS 

Several works on weak gravitational lensing reach 
up to the expressions that relate the optical scalars 
with the curvature components in terms of the tetrad 
adapted to the motion of the photons; we have here 
also presented expressions for the bending angle in 
terms of the curvature components. Furthermore, 
we have presented above expressions for the optical 
scalars and deflection angle directly in terms of the 
matter components of the sources of the gravitational 
lens, valid for an extended class of matter distribu- 
tions. Our expressions circumvent several deficien- 
cies as are: gauge dependence, lack of explicit expres- 
sions, neglect of spacelike components of the energy- 
momentum tensor, etc. The extension of this study 
to sources with different structure and to the cosmo- 
logical background will be presented elsewhere. 
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